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Abstract 

In these lectures we introduce the functional renormalization group out of equilibrium. While in thermal equilibrium 
typically a Euclidean formulation is adequate, nonequilibrium properties require real-time descriptions. For quantum 
systems specified by a given density matrix at initial time, a generating functional for real-time correlation functions 
can be written down using the Schwinger-Keldysh closed time path. This can be used to construct a nonequilibrium 
functional renormalization group along similar lines as for Euclidean field theories in thermal equilibrium. Impor- 
tant differences include the absence of a fluctuation-dissipation relation for general out-of -equilibrium situations. The 
nonequilibrium renormalization group takes on a particularly simple form at a fixed point, where the corresponding 
scale-invariant system becomes independent of the details of the initial density matrix. We discuss some basic exam- 
ples, for which we derive a hierarchy of fixed point solutions with increasing complexity from vacuum and thermal 
equilibrium to nonequilibrium. The latter solutions are then associated to the phenomenon of turbulence in quantum 
field theory. 

©2012 Published by Elsevier Ltd. 
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1. Introduction 

Thermal equilibrium properties of many -body or 
field theories are known to be efficiently classified 
in terms of renormalization group fixed points. A 
particularly powerful concept is the notion of in- 
frared fixed points which are characterized by uni- 
versality. These correspond to critical phenomena 
in thermal equilibrium, where the presence of a 
characteristic large correlation length leads to in- 
dependence of long-distance properties from de- 
tails of the underlying microscopic theory. In con- 
trast, a classification of properties of theories far 
from thermal equilibrium in terms of renormaliza- 
tion group fixed points is much less developed. 



The notion of universality or criticality far from 
equilibrium is to a large extent unexplored, in par- 
ticular, in relativistic quantum field theories. Here, 
the strong interest is mainly driven by theoretical 
and experimental advances in our understanding 
of early-universe cosmology as well as relativistic 
collision experiments of heavy nuclei in the labo- 
ratory. 

In the latter contexts, a particular class of non- 
thermal fixed points has attracted much interest in 
recent years. It is associated to the phenomenon 
of turbulence in quantum field theory, where a uni- 
versal power-law behavior describes the transport 
of conserved quantities. Traditionally, turbulence 
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is associated mostly with the dynamics of vortices 
in fluids HI] but also nonlinear waves can show tur- 
bulent behavior fl. In particular, it is well-known 
that interacting quantum field theories can lead to 
nonlinear dynamics and wave turbulence, even for 
very weakly coupled theories. 

Among the best studied theoretical examples in 
relativistic quantum field theory are scalar inflaton 
models for the dynamics of the early universe 
In a large class of models, the strongly accelerated 
expansion of the universe after the Big Bang is fol- 
lowed by turbulent behavior of the inflaton field be- 
fore thermal equilibrium is achieved. The connec- 
tion to the well-established phenomenon of weak 
wave turbulence in the presence of small nonlin- 
earities has been studied in great detail fl. Here, 
weak wave turbulence is associated to an energy 
cascade from small to high wave numbers. Only 
recently it has been realized that the direct energy 
transport towards higher wavenumbers is part of a 
dual cascade, in which also an inverse particle flux 
towards the infrared at small wave numbers occurs 
@|- One of the striking consequences is Bose con- 
densation far from equilibrium I6J. Similar scal- 
ing phenomena may also occur for gauge field dy- 
namics in the context of heavy-ion collisions at 
sufficiently high energies iM- llOll . or also for the 
nonrelativistic dynamics of ultracold atoms |[Tll - 
1311 . The emergence of same macroscopic scal- 
ing phenomena from very different underlying mi- 
croscopic physics is a formidable manifestation of 
universality far from equilibrium. 

Understanding the dominant collective phenom- 
ena in quantum field theories far from equilibrium 
represents a major challenge. Important phenom- 
ena, such as the infrared particle cascade and sub- 
sequent Bose condensation mentioned above, are 
genuinely nonperturbative and require suitable ap- 
proximation techniques. Here a nonequilibrium 
functional renormalization group approach 
2 in . or related real-time functional integral tech- 
niques based on n-particle irreducible (nPI) ef- 
fective actions B22I1 . can serve as a very useful 
means to gain analytic understanding. Other im- 
plementations of the renormalization group idea 
in this context include the so-called numerical 



renormalization group approach [23]], the time- 
dependent density matrix renormalization group 
1 2411 . the real-time renormalization group in Liou- 
ville space 12511 . or flow equations describing in- 
finitesimal unitary transformations [26]. These ap- 
proaches can be complemented by numerical simu- 



lations in (classical-statistical) nonequilibrium lat- 
tice theories 1 27 - 3^] or using kinetic descriptions 
113 111 in their respective range of applicability. 

In these lectures we discuss basic properties of 
the nonequilibrium renormalization group for the 
scale dependent generating functional of 1PI corre- 
lation functions in relativistic quantum field theory, 
following closely Ref. 11811 . Concentrating on the 
explicit example of a A^-component scalar field the- 
ory allows us to focus on the relevant differences to 
Euclidean treatments in vacuum or thermal equilib- 
rium, without introducing too much formalism and 
abstract notation. With the help of standard refer- 
ences from the functional renormalization group in 
Euclidean spacetime I32H42I1 . one can apply these 
considerations in a similar way to include fermions 
or gauge fields. 

We begin in section |2]by emphasizing some im- 
portant differences between thermal equilibrium 
and nonequilibrium. This will help us to under- 
stand why there exists a hierarchy of fixed point 
solutions with increasing complexity from vacuum 
and thermal equilibrium to nonequilibrium. In sec- 
tion [3] we introduce the notion of scaling behavior 
far from equilibrium in the context of weak wave 
turbulence for the description of stationary trans- 
port of conserved quantities. The nonequilibrium 
functional renormalization group is introduced in 
section |4j where we derive the relevant flow equa- 
tions. In section[5]we solve the flow equations in an 
approximation based on a resummed large-A^ ex- 
pansion to next-to-leading order. The results are 
used in section [6] to determine scaling exponents 
for nonthermal fixed points. We summarize and 
give an outlook in section [7] 



2. Thermal equilibrium vs. nonequilibrium 



2.1. Statistical and spectral functions 



All information about a quantum theory is en- 
coded in its correlation functions for a given den- 
sity matrix q. In thermal equilibrium, for the case 
of a canonical ensemble with inverse temperature 
/3 = l/T and Hamilton operator H, the density ma- 
trix is given by p (eq) ~ e~P H and it is normalized 
such that Tr g (eq> — 1 . A thermal real-time correla- 
tion function for a Heisenberg field operator <D(x) 
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is given by the time-ordered trace 

G^\x - y) = Tr [ e <*> T O(*)(D(v)} 

= <O(x)<&(y)> eq 0(* -/) 
+ (<!>(y)®(x)} eq 0(y Q - x°) , (1) 

for a two-point function, and involves n fields for 
a n-point function. Here, To is the time-ordering 
operator, x = (x°,x) denotes the time x° and space 
x variable^] and, to be specific, we consider real 
scalar fields. 

Nonequilibrium typically requires the specifica- 
tion of a density matrix g(to) at some initial time 
fo where g(to) + g^. The task of nonequilibrium 
quantum field theory is then to determine the real- 
time evolution of correlation functions such as the 
two-point function 

G(x,y) = Tr< £ .(f )T (D(x)<D(v)> = <T <D(x)O(v)> , 

(2) 

for times x°,y° > to- In general, translational in- 
variance does not hold out of equilibrium and thus, 
two-point functions will depend on both x and y 
separately, i.e. 

G(x, y) = F(x, y) - l -p(x, y) sgn(x° - /) . (3) 

In the second line of (O we introduced the sta- 
tistical two-point function F(x, y) and the spectral 
function p(x, y). These are given by the expectation 
values of the anti-commutator of the scalar fielc|^| 

F(x,y)=^(mx),<5>(y)}) , (4) 



and the commutator 



fKx,y) = i([<SKx),<SXy)]) , 



(5) 



respectively. The decomposition Q follows from 
(Q]i with the elementary properties of the Heaviside 
step function, 6(x°-y )+8(y°-x ) = 1 andsgn(jc°- 
y>) = 8(x° -y°) - 6(y° - x°). The spectral function 
is also related to the retarded and advanced Green's 
function by G R (x,y) = 0(x° -y°)p{x,y) = G A (y,x), 
respectively. Loosely speaking, the spectral func- 
tion p(x,y) determines which states are available 
while the statistical function F(x,y) contains the 
information about how often a state is occupied 
122]. 



1 We use a metric with signature (+, — , -, — ), 

2 In this section, we will assume for simplicity that we are in 
the symmetric phase where the field expectation value <p(x) = 
(<t>(X)) vanishes. Otherwise the connected statistical function is 
given by F(x, y) = \ <{$(*), *(y)(> - #x)0(y). 



2.2. Absence of a fluctuation-dissipation relation 

In contrast to the general nonequilibrium case, 
it is an important simplification of thermal or vac- 
uum theories that the statistical © and spectral 
function © are related by the so-called fluctuation- 
dissipation relation. We discuss this relation for 
the two-point correlation function with a canonical 
density matrix g (eq> ~ e^ H . For times t > the 
thermal two-point function ([TJ reads 

Tr{e-P H <!>(t,x)<l>(0,y)) 

= Tr [e~ pH e /iH (D(0,3')e^ H ®(t,x)\ . (6) 

In the second line we have used the invariance 
of the trace under cyclic changes and inserted 
e -P H eP H - i , From the real-time Heisenberg evo- 
lution <&(t,x) = e iH '<&(0,x)e~' H ', we may define 
in complete analogy <$(-ij3,x) = eP H <&(0,x)e~P H 
as an extension to 'imaginary times' l43ll . We 
can state the above relation directly in terms of 
the statistical (O and spectral function <(3j using 
<<D(f,x)O(0,j0> = F(t,0;x,y) - ip(t,0;x,y), and 
correspondingly for the r.h.s. of (|6). Since equi- 
librium is spacetime translation invariant, (|6) then 
reads 



F (eq) - y) - -p ieq) (x - y) 



y) 



/>=0 



~(eq) 



(7) 



Translation invariance also makes it convenient 
to consider the Fourier transform with real four- 
momentum p = (p°,p), that is 



(8) 



and equivalently for p (eq \x - y). Taking for a mo- 
ment for granted that these integrals can be prop- 
erly regularized and defined for the considered 
quantum field theory, we obtain from (0 

F (e(i) ( P ) - 1 -P {n \p) 

= e^ p (^(p) + |p<"»(p)) . 

This can be written in the form of a fluctuation- 
dissipation relation as 

F^\p) = -i (n BE (/) + i)p<"»(p) , (9) 
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where 

«be(/) = ^ , (10) 

is the Bose-Einstein distribution function. Since 
this distribution depends on frequency p° and 
not on spatial momenta p, it relates the anti- 
commutator expectation value of fields (F) and 
the respective commutator (p) in a nontrivial way. 
Moreover, at zero temperature the distribution 
function hbe is zero and the anti-commutator and 
commutator are directly proportional to each other. 
In contrast, for a nonequilibrium density matrix 
no such constraints exist in general. As a con- 
sequence, F and p are linearly independent out 
of equilibrium. In particular, the absence of a 
fluctuation-dissipation relation will allow us to ob- 
serve new types of scaling solutions, beyond those 
known from thermal equilibrium or the vacuum. 

2.3. Infrared fixed points 

Renormalization group fixed points correspond 
to scaling solutions for correlation functions. Be- 
fore computing them from first principles in the 
main part of these lectures, we illustrate here 
heuristically possible scaling behaviors in and 
out of equilibrium. For that purpose, we con- 
sider translationally invariant spectral and statis- 
tical two-point functions, F(x - y) and p(x - y), 
assuming also spatial isotropy to limit the number 
of possible scaling exponents. The Fourier trans- 
forms, F(p°,p) and p(p°,p), then depend on real 
frequency and momenta, where any scaling ansatz 
has to take into account the different possible scal- 
ings of spatial momenta vs. frequencies. This dif- 
ference is described by the 'dynamical critical ex- 
ponent' z. Furthermore, we denote the overall scal- 
ing of the statistical two-point function by the 'oc- 
cupation number exponent' k and of the spectral 
function by the 'anomalous dimension' rj. The 
scaling behavior may then be described as 

F(p°,p) = s 2+K F(s z p°, sp), (11a) 
p{p\p) = s 2 - T 'p(s z p\sp), (lib) 

for any real scaling parameter s > 0. 

Before we consider such a scaling behavior for 
the nonequilibrium case, it is very instructive to 
first insert (111a) and (II lbb into the fluctuation- 
dissipation relation (0 in order to see how this 
constraint relates the different exponents. We re- 
strict the discussion to frequencies (and momenta) 



much smaller than the temperature T. This would 
be the relevant range for scaling behavior, for in- 
stance, near second-order phase transitions in ther- 
mal equilibrium. In the infrared, where p° <K T, 
we see that the distribution function dTOb assumes 
the scaling form 

«be(p°) ~ ^ • ( 12 ) 
p" 

In particular, occupation numbers become large 
such that the 'proportionality factor' «be(p ) +1/2 
appearing in (0 can always be replaced by (TTZt for 
sufficiently small p°. Using dl lab and J 1 lbb we 
can then directly read off from <j9j the values of k 
for the vacuum and thermal cases: 

Vacuum (T = 0) : k = —tj , (13a) 
Thermal (T + 0) : k = -tj + z. (13b) 

As was mentioned in the introduction for the ex- 
ample of wave turbulence, also far from equilib- 
rium there exist important scaling solutions. These 
solutions can be spacetime translation invariant, 
however, they are in general not constrained by a 
fluctuation-dissipation relation. In that case, we 
may always write down a relation of the form 

F(p) = -i(n(p)+^p(p), (14) 

with some generalized 'distribution function' n(p) 
defined from the ratio of F(p) and p(p). However, 
in contrast to (fTOt for the case of thermal equilib- 
rium, here n(p) will in general depend both on fre- 
quency p° and spatial momentum p. Then from 
(11 1 at and (11 lbb itself no relation between the ex- 
ponents follows, but only that the distribution func- 
tion n(p) scales as 

n(p°,p) = s K+n n(s z p°,sp) , 

for n(p) » 1 /2. It will require an actual calcula- 
tion to determine a - at a nonthermal fixed point, and 
we will show from the nonequilibrium renormal- 
ization group in a large-A^ approximation to next- 
to-leading order that possible scaling solutions are 
SB 

Nonthermal: k — -tj + z + d , (15a) 
K = -j] + 2z + d. (15b) 

Here, d is the spatial dimension, and we will see 
that they describe the phenomenon of strong turbu- 
lence associated to particle (115al i and energy (115bb 
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cascades, respectively SUB 14411 . Because the 
dimensionality of space enters, the nonthermal val- 
ues for the exponent k and the corresponding fluc- 
tuations can be very large depending on d, which 
has been confirmed also in lattice field theory sim- 
ulations J5 , 13 , 44 , 45 ] . Comparing the nonthermal 
results with (II 3ab and d 1 3bb , one observes a hierar- 
chy of possible fixed point solutions with increas- 
ing complexity from vacuum, and thermal equilib- 
rium, to nonequilibrium. 

3. Basics of stationary transport 

3.1. Boltzmann transport 

Scaling behavior far from equilibrium is typi- 
cally discussed with the help of kinetic theory or a 
Boltzmann equation, which can describe transport 
properties of dilute, weakly interacting many-body 
systems. To make contact with the literature, we 
discuss the basic concepts of stationary transport in 
that way before we start from the nonequilibrium 
renormalization group in quantum field theory in 
section |U 

A Boltzmann equation describes the dynamics 
in terms of a single-particle distribution function 
n(t,p), which depends on time t and spatial mo- 
mentum/; for spatially homogeneous systems. The 
rate of change in the distribution of particles equals 
the difference between the rates at which particles 
in a phase space region are generated or lost due to 
collisions. For bosons there is an enhancement of 
the rate if the final state is already occupied. The 
collision terms may be expressed in terms of scat- 
tering cross sections and distribution functions. A 
Boltzmann equation describing 2 <-> 2 scatterings 
involving four particles reads 

d,n p (t) = c/r p i„23 \{n p + 1) («7 + \)n 2 n 3 

Jl,2,3 

- n p m (n 2 + 1) (n 3 + 1)] = C 2 ^i(t,p) , 

(16) 

where Jj 2J c/r p i„23 denotes the relativistically in- 
variant measure to be specified below. We have 
written n p = n(t,p) and «,■ = n(t,ki) as a compact 
notation for the distribution functions. 

Equation ( TToT ) can be understood as arising from 
the lowest-order perturbative contribution of the 
respective quantum field theory ]22|]. Since the 
functional renormalization group treatment start- 
ing with section|4]will include also the perturbative 



behavior, we only give here some relations to facil- 
itate comparisons with the literature. We consider 
the example of a relativistic, real scalar field theory 
with mass m and quartic self-interaction A, whose 
Lagrangian density is 



£ = l -(dM - -m 2 (D 2 - -CD 4 
2 M 2 4! 



(17) 



We introduce center and relative (Wigner) coordi- 
nates 

v x + y 

X — , r — x — y , 



and Fourier transform both the statistical (|4) and 
spectral function <j3j with respect to the relative co- 
ordinates 



F(X,p) = J", 
p(X,p) = J", 



d d+1 r e ipr F (X + r/2,X - r/2) , 
d d+1 re ipr p(X + r/2,X - r/2) . 



The Fourier transform of the statistical function is 
real and that of the spectral function is purely imag- 
inary due to their anti-commutator and commutator 
definitions, respectively. 

For spatially homogeneous systems, the correla- 
tion functions only depend on time t = X° and four- 
momentum p. We can define a time-dependent 
'distribution function' n(t,p), depending on four- 
momentum p, by writing 

F{t,p) = -i\n{t,p) + ^p{t,p), (18) 

which reduces to ( TPfl i for time translation invari- 
ant systems. If the spectral function is taken to be 
of the translation invariant lowest-order (free-field) 
form 

p (0 \p) = 2m sgn (p°) 6((p ) 2 - a? p ) , (19) 
with single-particle energy a> p , then 

X°° dD° 
-j^2p Q p (0 \p)n(t,p), (20) 

corresponds to the distribution function employed 
in the Boltzmann equation ( TToT ). The definition ( |20| ) 
ensures that the single-particle distribution func- 
tion n p (t) is evaluated for on-shell four-momentum 
with only positive energy, i.e. p° = a> p . The rel- 
ativistically invariant measure appearing in (Q~6} is 
then, for the theory dTTl . given by 

J w pl 6 2co„ J [\ = \{2n) d 2o H/ 

x (27T) d+,i 6 {d+Y \p +h-k 2 - h) , 

(21) 
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where to lowest order oj p = ~\Jp 2 + m 2 and we write 
a>i = a)(ki). 

More precisely, in a gradient expansion to low- 
est order in the number of derivatives with respect 
to the center coordinate X° and in powers of the 
relative coordinate r°, the spectral function for spa- 
tially homogeneous systems obeys 1 22] 



(22) 



in agreement with the constant free-field form (TT9V 
In contrast, the statistical function F{t,p) to this 
order can be time-dependent and its evolution is 
given by 



I 

Jo 



dn° 

JL2p°d l F(t,p) 
2tt 



X°° do 
■j^2p°p p d t n(t,p) = C(t,p), (23) 



which defines the collision term C(t,p) describing 
the effects of interactions to lowest order in the gra- 
dient expansion. The restriction of the collision 
term to 2 «-> 2 scatterings, which are the leading 
perturbative contributions ~ A 2 , gives the Boltz- 
mann equation ( TTol l. 

We emphasize that the Boltzmann equation has 
a limited range of applicability. In particular, it 
is restricted to the perturbative regime. Even for 
weak coupling, A <k 1, nonperturbative correc- 
tions can play a crucial role once the momentum 
modes become highly occupied. More precisely, 
equation ( TTol ) can only be applied for parametri- 
cally small occupancies, n p <k l/A. Nonpertur- 
batively large occupation numbers n p ~ l/A lead 
to important corrections, since the 'thin-gas' ap- 
proximation underlying the Boltzmann equation no 
longer holds and multiparticle scatterings have to 
be taken into account. The calculation of these 
nonperturbative corrections will be a central topic 
for the nonequilibrium functional renormalization 
group below. However, for the moment we will 
consider what one expects from perturbation the- 
ory. 

More precisely, we will analyze (fTol i in 
the 'classical' regime of occupation numbers 
l/A » n p » 1, where the collision term can be ap- 
proximated by 

C { ^ 2 (t,p) = c/r p i„23 \(n p + ni)n 2 n 3 

Jl,2,3 

- ripii!^ + «})] . (24) 



inverse particle 
cascade to the IR 



n(p)~l/A 



direct energy cascade 
to the UV 



l/A > n(p) > 1 



quantum/dissipative 
regime 



\ 

X 


















log p 



Fig. 1 . Illustration of the dual cascade for scalar quantum field 
theory. Different scaling exponents can occur in different mo- 
mentum regimes, which are associated to stationary transport of 
conserved quantities. 



For n p < 1 quantum or 'dissipative' processes will 
play an important role, which will obstruct scal- 
ing at sufficiently high momenta. The situation for 
the real scalar field theory is schematically sum- 
marized in Fig. [1] which sketches the occupation 
number as a function of momentum on a double- 
logarithmic scale such that straight lines corre- 
spond to power-law behavior. Different slopes can 
occur in different momentum regions. The fact that 
each scaling exponent is associated to an approxi- 
mately conserved quantity, such as energy or par- 
ticle number, will be explained next for the pertur- 
bative regime using the Boltzmann equation. 

3.2. Weak wave turbulence 

We may write the Boltzmann equation (TTol l for- 
mally as a continuity equation 



dt£ + V p • j p = , 



(25) 



for the energy density e(p, t) = to p n(p, t) in mo- 
mentum space, where the divergence of the energy 
flux j p is given for the collision term (l24l by 



Since total energy is conserved for the relativistic 
quantum field theory, we may ask for its effect on 
the dynamics. We consider the case of an isotropic 
system where the only nonvanishing part of the 
flux is given by its radial component. Integrating 
the continuity equation d25l ) over the volume of the 
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sphere B(k) of radius k in momentum space, we 
obtain 



JB(k) 



d d pd,e = -(2n) d A(k) 



(26) 



which states that the change of net energy con- 
tained in B(k) is given by the flux A(k) through the 
boundary dB(k): 



-(2n) d A(k) 



2n d ' 2 

r(d/2) 



f 

Jo 



d\p\\p\ d -'u p CfUt,p) 



(27) 



Thermal equilibrium is characterized by the van- 
ishing of the collision term and a zero net flux 
A(k). Here, however, we are interested in possi- 
ble stationary solutions of ([Tol l where the distribu- 
tion n p characterizes some steady state that is far 
from equilibrium, known as weak wave turbulence 
J3. For such a steady state to exist, the distri- 
bution function n p needs to satisfy the stationarity 
condition 

c£> 2 [n] = . (28) 

However, in contrast to thermal equilibrium, it 
has a nonvanishing flux. Such a 'flux state' de- 
scribes the stationary transport of conserved quan- 
tities 10]. 

Here, we show that these stationary states corre- 
spond to the situation where A(k) becomes scale- 
independent. To investigate the behavior of the en- 
ergy flux under scaling transformations, we make 
for the single-particle energy a scaling ansatz 

u>(sp) = su>(p) , 

using the linear dispersion for the relativistic scalar 
field theory ( fT7T > at sufficiently high momenta 
p 2 » m 2 . Similarly, the occupation number distri- 
bution is taken to obey the scaling form 

n(sp) = s~ K n(p) , 

with the occupation number scaling exponent k in- 
troduced already in section [2] Together with the 
scaling property of the measure, 



J ■ 

J 1,2,3 



dr p i„23(sp, ski, sk 2 , ski) 

! 

= f dr pU + 23 (p,kuk 2 M) , (29) 

Jl,2,3 

with fx 4 = (3d - 4) - (d + 1) = 2d - 5 for the scalar 
field theory case d2Tl i, we see that the collision in- 
tegral satisfies 



This result can easily be generalized for m-particle 
scattering processes 

where C,„ denotes the collision term for the case 
of m-particle scattering in the classical regime, and 
yu m = {m - 2)d -m — l. Using, accordingly, u>(p) = 
\p\co(l) and n(p) = \p\~ K n{\) etc. for the isotropic 
system, the flux (l27T i can be seen to give 



A(k) 



w(l)Cf(l) 



d + 1 — K(m — 1) + fi„ 



c l^d+\-K(m-l)+fi„, 



(30) 

For A(k) to become independent of the scale k up 
to logarithmic corrections, the ^-exponent d + 1 - 
K(m— l)+fi m must vanish. For the scaling exponent 
k this gives 

Mm +d+l m 

k= — =d -. (31) 

m - 1 m - 1 

Since the denominator of (f30b also vanishes in this 
case, the limit 



C 



(cl) 



->d-m/(m-l) d + 1 — K(m - 1) + jl n 



- const. 



has to exist. Therefore, the collision integral must 
have a zero of first degree in d+ 1 -K(m- l)+j-i m 118 , 



4611 . Here, we will assume that the stationary state 
exists and refer to the literature for further details. 
Specifically, for the scalar theory with quartic self- 
interaction in d — 3 spatial dimensions, we have 
the scaling exponent 



d = 3 



for the transport of energy over some range of mo- 
mentum scales. This is the so-called energy cas- 
cade, which is well-known from the perturb ative 
theory of weak wave turbulence flS]- We also 
note that from the relativistic scaling ansatz one 
observes that the scaling exponent associated to 
momentum conservation is the same as for the en- 
ergy cascade. 

If the dynamics is approximated by the Boltz- 
mann equation ( TTBT l. i.e. if only the perturbatively 
leading 2 <-> 2 scattering processes are taken into 
account, then particle number is conserved. Of 
course, there are total particle number changing 
processes in the considered relativistic quantum 
field theory. However, these processes appear at 



s 



/. Berges et aL / Nuclear Physics B Proceedings Supplement 00 (2013) 1 424\ 



higher order in the coupling, such that for A <<c 1 
inelastic scattering rates are much smaller than the 
elastic ones 14711 . Because of this separation of 
scales, there can be important consequences of 
approximate particle number conservation for the 
phenomenon of weak wave turbulence. Since the 
perturbative particle flux is simply given by the 
momentum integral of the collision integral C 2 ^ 9 , 
going through the corresponding steps as above 
one may easily verify that the scaling exponent 



k = d ■ 



m + 1 
m — 1 



describes the stationary transport of particles. For 
the special case of m = 4, and d — 3, this particle 
cascade is characterized by the exponent 



d = 3 



K — — 

3 



So far, we have only considered effects of 2 — » 2 
particle scattering. For stationary turbulence in 
quantum field theories, however, the dynamics can 
lead to a nonvanishing field expectation value or 
Bose condensation far from equilibrium fl. This 
gives rise to an effective three-vertex (m = 3) such 
that the energy cascade in three dimensions is char- 
acterized by the scaling exponent J3l 



d = 3 



3 

K = — 

2 



In these lectures, we will consider only the case of 
a vanishing field expectation value for simplicity. 

From the above discussion we have observed 
that there are two types of perturbative power-law 
distributions corresponding to stationary transport 
of energy and particle number, respectively. These 
may be realized in different regions of momentum 
space. However, the simple analysis cannot de- 
termine in which regions of momentum space the 
different scaling solutions are realized, or whether 
the cascades describe transport from small to large 
wavenumbers or vice versa. A thorough discus- 
sion leads to the picture of a dual cascade in scalar 
quantum field theory as illustrated in Fig. Q] 0, 
I2l LL2l |45[] • Its quantitative description at low mo- 
menta requires, however, to go beyond perturba- 
tion theory. It has been found that the approximate 
conservation of an effective particle number can 
be applied to the nonperturbative regime of scalar 
field theories at low momenta, which leads to dif- 
ferent values of turbulent scaling exponents than 



the perturbative analysis suggests QUil]. Thus, we 
are in need of reliable nonperturbative techniques 
to access all characteristic momentum regions of 
the nonequilibrium dynamics. In the following, we 
will consider the functional renormalization group 
as an ideal tool to access a wide range of scales in 
a unified framework. 

4. Functional renormalization group 

We want to obtain information about scaling so- 
lutions for correlation functions far from equilib- 
rium. In the previous section, starting from pertur- 
bative kinetic theory or the Boltzmann equation, 
we illustrated how such a scaling behavior may 
arise in a weakly coupled scalar theory. However, 
in the infrared, where occupation numbers become 
large and the system is strongly correlated, the per- 
turbative description is insufficient. Thus, we need 
a general framework to calculate correlation func- 
tions from first principles. Here, we write down 
a generating functional for nonequilibrium corre- 
lation functions 1221 14811 . which serves as a start- 
ing point to define the nonequilibrium functional 
renormalization group on a closed time-path fol- 
lowing the presentation of Ref. 118 1. 



4.1. Generating functional 

All information about a nonequilibrium quan- 
tum field theory can be efficiently described in 
terms of the nonequilibrium gene rating functional 
for correlation functions [ 2^ 48 ] . For given den- 
sity matrix go = g{to) at some initial time to, corre- 
lation functions can be obtained from the generat- 
ing functional 



Jx,C 



Z[J,R;q ]=Tiq T c exp/ { \ 

®(x)R(x,y)®(y)} , 



2 X.v.c 



for a scalar field theory as described by (TPTt in 
the presence of sources J and R. The introduction 
of the bilinear source term ~ R will be convenient 
for the derivation of the functional renormalization 
group equation, which is explained below. Here the 
time integration is taken over a closed time path C, 
i.e. l c = j c dx° f d d x, displayed in Fig. |2] SHI- 
The closed time-path appears because we want to 
compute correlation functions, which are given as 
the trace over the density matrix with time-ordered 
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to 



c + 



1^ 



Fig. 2. Closed time path C. 



products of Heisenberg field operators, as exempli- 
fied in section |2 Representing the trace as a path 
integral will require a time path where the initial 
and final times are identified, which is discussed in 
more detail below. 

Above, Tc denotes time-ordering along the con- 
tour C. As seen from Fig. [2] this contour con- 
sists of an upper C + and a lower branch C~ where 
the time-ordering on the lower branch is reversed. 
To extract correlation functions efficiently, the field 
$>(x) may be written in terms of <l> ± (x ,jt:) where 
the ±-index denotes on which part of the contour 
C ± the time argument is located. E.g. the contour 
integration for the source term in (l3Zt takes the 
form 



I 



<&(x)J{x) 



= f (D + (x)7 + (x)+ f (D -(x)J-(x) 

Jx,C* Jx.C- 

=1 



d d x{<S> + (x)J + (x)-<S>-{x)J-{xj) , 



where the minus sign comes from the reversed 
time-ordering along C~. Setting the sources J and 
R to zero in (l32l i we obtain the partition sum 



Z[J,R;qo]\ 



Trpo = l 



from the normalization of the density matrix. 

For instance, taking the second functional 
derivative of the generating functional (l3Zt with re- 
spect to the classical source J + defined on C + and 
setting all sources J and R to zero, we obtain 



5 2 Z[J, R;q ] 



i6J + {x) i5J + (y) 



= <T O(*)O(y)) = G ++ (x,y). 



j,r=o 



Here we have used that Tc is identical to standard 
time ordering To in this case, since both x° and 
y° lie on the upper part of the contour, i.e. on C + . 
We also introduced the notation G ++ {x,y) in order 
to distinguish this correlator from the other pos- 
sible second functional derivatives with respect to 
the sources J + ,J~ and setting J, R = afterwards. 



These can be written as: 

G ++ (x,y)^(<b(x)<b(y)0(x°-y o ) 

+ <DCv)(D(x) 0(y° - x )) , 
G-(x,y) = (®(x)<l>(y)e(y - x°) 

+ ®(y)®(x)9(x -y°j) , 
G + -{x,y) = mymx)) , 
G- + {x,y) = (cD(x)<B(y)> . 

We emphasize that not all of the above two-point 
functions are independent. In particular, using the 
property 0(x° -y°) + 6(y° - x°) = 1 of the Heaviside 
step function one obtains the algebraic identity: 

G ++ (x, y) + G—(x,y) = G + -(x,y) + G- + (x,y) . 

(33) 

This identity will be of use later on. 

4.2. Functional integral 

To simplify the evaluation of correlation func- 
tions we write the generating functional (l32t in 
terms of a functional integral representation. For 
an intuitive presentation we follow standard tech- 
niques (see e.g. lf5ll0 : We evaluate the trace using 
eigenstates of the Heisenberg field operators cD* at 
initial time fo, 

®*(t ,x)\<p ± ) = i p t (x)\<p ± ), 

such that (l32l l may be written as 

Z[J,R;q ] 

= j [d<p + ] {<p + I go T c exp i j jT ®(x)J(x) 
+ \ f <S}{x)R{x,y)®(y)\\ip + ) . (34) 

1 Jx,y.C ) 

Here the integration measure is given by 

j [dipt] = J IJd^ix). (35) 

With the insertion 

[d< Po -]\<p-)(<p-\ = l, (36) 
we may bring ( f34b to a form 



Z[J,R;g ] = j [# ( ;][# ] (<P + \go\<p-) 

x (<p-,t \<p\t Q ) JR . (37) 
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Here the transition amplitude in the presence of the 
sources is given by 



(ip ,t \<f + ,t ) 



J.R 



= (<f-\T c «p» jjT <D(. 

+ 5 X c ® {x)R(x ' y)<S>{y) } 1 * +) ■ (38) 

This matrix element can be written as a functional 
integral over the fields ip ± 

(ip-,to\ip + ,t ) JR 

= j [d<f + ]'[d(f-]' espils[<p} + f VteW*) 

+ ^X c ^ x)R{x > y)ip(y) } ' (39) 



which is essentially the same procedure as em- 
ployed to obtain standard path integral expressions 
for vacuum of equilibrium matrix elements 1143 1. 
Here S [if] is the classical action for the scalar the- 
ory and the measure is given by 



(40) 



x; x u > to 



Here, the functional integration goes over the field 
configurations ^ ± (x°,jr) that satisfy the boundary 
condition ^(x = to,x) = (f^ (x). 

The above expression d37l > together with 
$3% displays two important ingredients entering 
nonequilibrium quantum field theory: the quan- 
tum fluctuations described by the functional inte- 
gral with action S that determines the transition 
matrix element, and the statistical fluctuations en- 
coded in the averaging procedure over the initial 
conditions as specified by the initial density ma- 
trix go. 

4.3. N-component scalar field theory 

So far, we have not specified any internal field 
degrees of freedom. In the following, we consider 
a TV-component vector field <t> a with a — 1, . . .,N 
for an <9(A0-symmetric theory. The number of field 
components will later serve as an expansion pa- 
rameter that allows for a controlled approximation 
of renormalization group equations. The classi- 
cal action for the 0(N)-model with a quartic self- 



interaction is given by 

S[<f] = - f if a (x)iD a l(x,y)if b (y) 

1 Jx.y.C 

~ J <Pa(x)<fa(x)iPb(x)<fb(x) , 

(41) 

where the time integration runs over the contour 
and iD~} denotes the free inverse propagator satis- 
fying 

(-□* - m 2 } iD ab (x,y) = 5 ab 6^ +l \x - y) . 

Using again the +-index to denote on which part 
of the contour C ± the time argument is located, the 
free part of the action takes the form 

Sol<P + ,<f-l = \^ (if + a (x)iD- a l(x,y)ift(y) 

- (f-(x)iD-^(x,y)ifl(y)) . 
The interaction term is written as 

S inl [ip + ,if~] = J (<ft(x)fa(x)<ft(x)<ft(x) 

~ <Pa(x)<Pa(x)<Pl(x)(fl(x)) . 



The time integration is implicitly defined along 
the positive branch of the contour, i.e. J = 

J t °° dx° J d d x and the minus sign in front of the ip~- 
components of the action comes from the reversed 
time ordering along the lower branch C~. 

4.4. Quantum vs. classical dynamics 

In order to discuss the different origins of quan- 
tum and of classical-statistical fluctuations, it is 
convenient to intraductal 

fa = J (<Pa + <Pa) > <P<> = Pa~ ~ fa ■ (42) 

For the rest of the lectures we will use this basis for 
our calculations. The action (l4"TT i takes the follow- 
ing form 

S [<f, ij>\ = S [if, y]+S int [if, if] , (43) 
where the free part is given by 



So\<P, <f>] 



J Va( 



(x)iD- a l(x,y)ip h (y) , (44) 



In order to prevent a proliferation of symbols, the notation 
does not distinguish the new definition of ip a from its earlier use 
since no confusion for the following can occur. 
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and the interaction part reads 

SmtW,<P\ = J ipa(x)<Pa(x)<Pb(x)<Pb(x) 



(x)(p b (x)<p b (x) , (45) 



Finally, the linear source term can be written in 
this basis takes the form 



J(cp + a (x)J:(x)-if-(x)j;(x)) 

= J (<Pa(x)J a (x) + ip a (x)J a (x)) , 
whereas the bilinear source term is written as 



1 



(*£(*),¥>-(*)) 



R++{x,y) -R + ab (x,y) 



l ab 

-R:Z(x,y) 



■/< 

U x,y 



(ifi a (x) , ifi a (x)) 



<£(y> 



We have added the indices R, A, F, P to the bilinear 
sources that suggest a relation to the retarded, ad- 
vanced, and statistical components in the new ba- 
sis. This connection will be made more explicit in 
the following subsections. 

The two types of vertices appearing in the in- 
teraction part (l45l l are illustrated in Fig. [3] To un- 
derstand their role for the dynamics, it is impor- 
tant to note that one can also write down a func- 
tional integral for the corresponding nonequilib- 
rium classical-statistical field theory. The stan- 
dard derivation of the latter employs that the clas- 
sical field equation of motion can be obtained from 
S [ip, ip] by functional differentiation with respect 
to ip, 



6S [<p, ip] 



5(pa{x) 



■ (O x + m 
A 



3N 



) <Pa(x) 

<p a (x)(pb(x)<pb(x) = , (46) 



where the derivative had to be evaluated for tp = 
to eliminate terms originating from the part of the 
interaction term (|43T > that is cubic in tp. This clas- 
sical dynamics can then be implemented as a func- 
tional integral using the representation of the 6- 
functional 



6S [<p, ip] I 



6ip a (x) \g>=o 



9 



<P 



ip (f 



Fig. 3. Classical (left) and quantum vertex (right) in the scalar 
field theory. 



with the help of the 'auxiliary' field <p and further 
steps involving the integration with respect to <p 



Il52l - l54ll . We emphasize that in the exponent on the 
r.h.s. of (|47| | only terms linear in ip appear. In par- 
ticular, the interaction term ~ ip 3 , appearing with 
(PES! in the functional integral of the quantum the- 
ory, does not occur for the classical theory. In 
summary, the generating functionals for correla- 
tion functions are very similar in the quantum and 
the classical statistical theory. A crucial difference 
is that the quantum theory is characterized by an 
additional vertex. In a diagrammatic language, if 
loop corrections involving only the classical ver- 
tex dominate over those involving the quantum ver- 
tex or a combination of both, then the quantum 
dynamics can be approximately described by the 
classical-statistical field theory. This has been an- 
alyzed in great detail in recent years for nonequi- 
librium phenomena J5, 55, 5(|, and we will come 
back to this point in the context of turbulence be- 
low. 



4.5. Connected one- and two-point functions 



The generating functional for connected correla- 
tion functions is given by 



W = -/InZ . 



(48) 



We may calculate field expectation values by func- 
tional differentiation with respect to the sources J 
and J, which we denote as 



6W - s 5W 

(p a (x) , — = 4> a (x) 



SJ a (x) 



6J a {x) 



for fixed J and R. The connected two-point corre- 
lation functions are given by the second functional 
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derivatives 



6 2 W 



SJ a (x)SJ b (y) 

5 2 W 
8J a (x)dJ b (y) 

6 2 W 
5J a (xW~ b (y) 

6 2 W 
6J a (x)6J b (y) 



= G A b (x,y), 
= iF ah (x,y) , 
= iF ab (x,y) , 



(49) 



where G RA are the retarded/advanced propagators, 
F is the statistical propagator, and F is the 'anoma- 
lous' propagator. We note that the retarded and ad- 
vanced propagators satisfy the symmetry property 

G A b (x,y) = Gf a (y,x), (50) 

and for the statistical propagators, we have 

F a b(x, y) = F ba (y, x) , F ab (x, y) = F ba (y, x) . 

These properties follow directly from the definition 
of the propagators in terms of the second functional 
derivatives with respect to / and /. The spectral 
function p is given by the difference of the retarded 
and advanced propagators 

Pab(x,y) = G R & (x,y) - G A A (x,y) , 

and G R b (x, y) = p ab (x,y)8(x°-y°). It is important to 
note that the anomalous propagator F vanishes in 
the limit where the external sources are set to zero. 
This is a consequence of the algebraic identity (l33l l. 
since F = G ++ + G - G + - G + ~, as one may 
readily check by changing basis. More generally, 
in the absence of sources, we have [51] 



6W I 
rV./,,(\ )•/./.« 



= Mx) = . 



6 2 W 



5J a (x)6J b {y)\Jj,R RAF - t =o 



iF ab (x,y) = , 



and, correspondingly, arbitrary functional deriva- 
tives of the generating functional with respect to J 
vanish in the absence of sources. 



4.6. Functional renormalization group 

A most convenient derivation of the functional 
renormalization group equation starts from the 
two-particle irreducible (2PI) effective action 15711 . 
The latter is obtained as a Legendre transform of 
the generating functional (f48b with respect to the 



linear and bilinear source terms J, R. One obtains a 
functional of the field expectation values 0, 0, and 
the propagators G R , G A , F, F: 

r 2P1 [cf>,0,G R ,G A ,FF] 

= W- jj<P a (x)Ja(x) + ~4> a {x)J a {x)) 

- \ f {Ratify) [Ux)Uy) - iG* b (x, y )] 

+ R* b (x,y) [j> a (x)(f> b (y) - iG^{x,y)] 
+ R^ b (x,y) [cf>a(x)cf> b (y) + F ab (x,y)] 
+ R F ab {x,y) [Ux)My) + Fab(x,y)] } , 

where the sources depend on the field expec- 
tation values and the propagators, i.e. / = 
/ [<p, 0, G R , G A , F, F] etc. 

We may partially undo the Legendre transform 
for the bilinear sources S R * A,f ' f to obtain an expres- 
sion for the one-particle irreducible (1PI) effective 
action in the presence of these sources, i.e. 



T [<p, j>, fl R - A = r 2PI - Tr [g r r r + G A R P 



+ iFR F + iFR F] 



(51) 



Here the propagators depend on the fields and 
bilinear sources, i.e. G R = G R [<f>, </>, R R A FF ] 
etc. The above equation is our starting point for 
the construction of the functional renormalization 
group. We take the bilinear sources to depend on 
some characteristic momentum scale k > 0, which 
renders the effective action scale-dependent, 

r t [0,«j] = r [m,s** fJI ]. 

Taking now the derivative with respect to the scale 
k, denoting 



T k [4>,4>] = k 



d T k [(/>,4>] 
dk 



we get from (|5TJ the renormalization group equa- 
tion for the scale-dependent effective action: 

f* = ~\ Tr {GfRf + GX + iF k R F + iF k R F k ) , 

(52) 

Here we have used the fact that Tipi in ( Bil l is inde- 
pendent of the bilinear sources Rf KF - F . The flow 



equation (1521 corresponds to the Wetterich equa- 
tion for the effective average action iBitl . however, 
now evaluated on the closed time path. 
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4.7. Cutoff functions 

With suitable cutoff functions, the effective av- 
erage action may be viewed as a coarse grained 
effective action, which includes all quantum- 
statistical fluctuations with characteristic momenta 
above the scale k. If the infrared cutoff scale k is 
sent to zero, one recovers the standard 1PI effective 
action with all fluctuations included, i.e. 

T k=0 [0, 4>]=T [<p, ~d>, ^ F ' p = 0] . (53) 

Therefore, the sources R^' ' F ' F must vanish in the 
infrared limit. With this property, the renormaliza- 
tion group equation may be used to describe the 
flow starting from some high momentum scale A, 
where the microscopic physics is characterized by 
some classical action S , i.e. 



Km 1*0,0 -5 [0,0] . 



(54) 



The renormalization group flow then interpolates 
between the classical action and the full quantum 
effective action that appears when all fluctuations 
have been taken into account. 

We have seen that the renormalization group on 
the closed time path requires the specification of 
the various cutoff functions R k ' ' ' , which seems 
to allow for a wider class of possible choices than 
in the corresponding Euclidean field theories. In 
order to discuss what constraints can be obtained 
from the requirement (l54l . it is useful to start 
from the - now ^-dependent - generating func- 
tional 



Z k [J, 



[dip] [dip] exp i 



[s [<p, &] 



ip a (x)J a (x) + ip a (x)J a (x)} 

j (<p a (x)R£ ab (x,y)<p b (y) 

+ <Pa(x)Rf ah {x,y)tp b (y) 
+ <Pa(x)R F ah {x,y)(p b (y) 

+ <p a (x)R F ab (x,y)(p b (y)}\. (55) 



Here we have hidden the averaging over the initial 
density matrix qq in the notation, since we will not 
be concerned with the special choice of the initial 
conditions for the present purposes. Later on we 
will comment on the role of the initial conditions. 



The affective average action (IBTT l can then be writ- 
ten as 

T A .[0,0] = -zlnZ* - J^cf> a (x)J a (x) + (p a {x)J a {x)\ 



4I> 



(x)R£ ab (x,y)fo(y) 



+ 4>a{x)RlJx,y)(j> b {y) 
+ (pa(x)R F kab (x,y)d> b {y) 

+ Ux)Rl ab (x,y)~My))- (56) 
In the following, we employ the notation 



6T k 



P fx) EE p ( x ) = 

k ' aV Scp a (x) ' k ' aV 5Ux) 
By functional differentiation of (l56l l with respect 
to the fields <p a and (j> a , we obtain the equations of 
motion for the fields <t> and d>: 



ri >fl (*) = -•/«(*) - [Rl ab {x,y)cp b (y) 



+ \Ria b (x,y)$ b (y) + l -4> b (y)Rl ba (y,x) 



(57a) 



K. ab (^y)h(y) 



+ 2^t*(*,y¥*(y) + ^Uy)Ri ba (y,x) 



(57b) 



With these we may eliminate the sources J a {x) and 
JJjx) from ([56l. expressing them in terms of the 
first functional derivatives of T k . We then shift the 
fields in the generating functional ( l55l l by 



tp , p — » <p + ip , 



observing that the measure is invariant under these 
transformations. That way, from ((56), we fi- 
nally arrive at the following functional integro- 
differential equation for the effective action: 

= S [(/>, 4>] - i In J [dip] [dip] exp i js [0 + tp, !f> + ip] 

- S [0, 0] + f{<Pa(x)Tl a (x) + <p a (x)lf a (x)} 

~ 2 L [ ifia< > x ^ R t,b( x >y')'Pb(y) 

+ >f a {x)Rf ab (x,y)(p b (y) 
+ (p a {x)R F ab (x,y)ip b (y) 

+ (p a (x)R F ab (x,y)y b (y)}\. (58) 
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With the following representation for the 6- 
functional 

<%] = JWlexpjz J <p a (x)<p a (x) 

and equivalently for S[ip], one observes that the 
property (|54| | can be efficiently achieved with a 
class of cutoff functions chosen as 



Rf£(x,y) = R k (-□,) 6(x - y)5 ab , 



k,ab 
k,ab 

with the property 



lim Ri 



(59) 



(60) 



This ensures that the bilinear source terms in 
act as (5-constraints in the limit k — > A, thus sup- 
pressing fluctuations in the fields <p and <p. More 
generally, possible nonzero R F k ' F should not be cho- 
sen to grow as fast as R^' A for k — > A in order to 
comply with (O OH- 

The choice of vanishing R FF greatly simplifies 
the structure of the flow equations. In this case the 
exact flow equation for the effective average action 
d52i i becomes 

f k [cf>, # = Tr [GfRf + GX) . (61) 

Furthermore, from the symmetry properties (l50l 
of the propagators Gf' A , and the cutoff func- 
tions Rf A , one observes that Tr {G?[0, fakf) = 



TrjGj [(/>, 4>]R A j holds. However, it proves to be 
convenient to keep both the retarded and advanced 
functions in doTT ) as it simplifies the diagrammatic 
rules that will be introduced below. 

4.8. Propagators 

Above we have derived the exact flow equation 
for the effective average action F k that depends on 
the retarded and advanced propagators Gf and G A . 
It remains to relate the two-point functions to func- 
tional derivatives of T k . Here, we want to illustrate 
how to obtain these relations starting from simple 
identities for the sources J and /. As an example, 
we may consider the following identity 



6 (d+l \x-y)6 ab 
= 6J a (x) C(6Ja< 

sJ b (y) J z \s<pc' 



,,{x) 6<t> c (z) 5J a {x) 4> c (z)\ 



c(z) SJ b (y) 



Uz) Jb(y) J ' 
(62) 



where we have used that J - J[(f>, 0] is a functional 
of the field expectation values. Using the equation 
of motion 



S*- = -J a (x) - f <„,,(x, 



y)4>biy) , 



we may write the functional derivatives 
6J a (x)/6cf> c (z) and 6J a (x)/6<f> c (z) in terms of 
second functional derivatives of the effective aver- 
age action. Here we have exploited the symmetry 
property of the cutoff functions Rf' A and the van- 
ishing of R FF to write (I57bb in a somewhat simpler 
form. Furthermore, the functional derivatives of 
the field expectation values 



ail) 



s 2 w k 



54> c (z) 



S 2 W k 



SMy) 6J c (z)6J b (y) 5J b iy) 5J c {z)5J b (y) 

can be expressed in terms of the fc-dependent gen- 
erating functional W k . That way, we may rewrite 

as 



6 (d+1 \x-y)6 ab 

■X{( " r " 



5 2 W k 



d(f> a (x)Scp c (z) K ' uc * "1 6J c (z)6J b (y) 

2 n s 2 w k ) 



6$ a (x)6<f> c (z)6J c (z)6J b (y)} 



(63) 



With d49l this can be written in a compact matrix 
form, 

(lf+Rf)Gf + lfiF k = -l, 
where we have used the notation 

pfr ( \ = 5 2 r k [<f>,<i>] 

«^y>- 5 ~ Ux)5Uy) ■ 

for the functional derivatives of the effective aver- 
age action. 

Starting from the remaining three identities for 
the functional derivatives of the sources, that is 



and 



sJbiy) 

5J a (x) 



= S {d+l) (x-y)6 ab , 



= 0, ^4=0, 



5J b (y) 6J b (y) 
we obtain the set of equations 



(lf+R A )G A +lfiF k = -t 
(if* +R*) iF k + if G A = 0, 
(lf+R A )iF k+ lfGf = 0. 
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Together with our result from above, this linear 
system can be solved for the propagators Gf ' A , F^, 
and F k : 



iF k = - 
iF k = - 



(if + Rf) - if (if + r a )~ if 

(if + r a ) - if (if + Rff 1 rf 

rf -(if +**)(rf) _1 (rf 
if _ (if + i«) (if J- 1 (if + ^) 



(64) 



where the propagators depend on the field expecta- 
tion values (p a and (f> a , i.e. Gf' A = Gf' A [(p,(f>], etc. 
We emphasize that the above is valid for the choice 
R FF = 0, and for nonvanishing R FF the propaga- 
tors take a different form IU8I0 . 



4.9. Diagrammatic s 

From the flow equation for the effective average 
action ( loTb . we can construct the flow equations for 
arbitrary n-point functions by functional differenti- 
ation. As an example, we consider 



I 6 2 Gf [</>,$] 6R 



l a (x) 5<t> b (y) 



-R 



+ —z- 



5 2 G A [cp, 



-R'r 



d(f) a (x) 6(f) h {y) 



involving functional derivatives of the retarded and 
advanced propagators. For the retarded propagator 
we have, e.g. 



6 2 G R 



k,ab 



[0,0] 



U k,ae L k,efcd U k,fb 



_ /->R t4<M»I> I p<W , pA^ -pH s^R 
U k,ac l k.efcd \ l k + K k ) fg 1 k,gh U k, 

+ U k,ae l k,ef \ l k k ) fg k,ghcd 

1 

k.lu 



hb 



X 



Vk +K k) kM ikJjGk.jb 

_ U k,a e l k,ef\ l k +K k) fg l k,ghcd U k,hb 

+ plus terms involving three-vertices 

This seems rather complicated at first sight. How- 
ever, after taking derivatives, we can set the sources 
J — J = to zero, which corresponds to evaluat- 
ing all expressions at the field configuration {<p = 
(f>o(k), = 0) that extremize the effective average 



action. Here <po(k = 0) can be nonzero in the case 
of spontaneous symmetry breaking. According to 
the equation of motions, this configuration fulfills 



with (I57ab for / = J = (and R F ' F = 0) 



(x) = . 



More generally, all the functional derivatives of the 
effective average action with respect to (p vanish at 
the extremum 115111 . i.e. 



&T k [(/>,4>] 



6(p{xi) 6(f>{x2) ■ ■ ■ S(p(x n ) 



= 



Mk),4>=0 



In particular, we have T k [4> = (po(k), <p — 0] = 
so that the d64l i simplify considerably. The anoma- 
lous statistical propagator F k [(p = (po(k),^> = 0] is 
zero and the nonvanishing propagators are given by 



Gtb = -(rt* + R A ) 



ah 



iF k , ab = (Gflf Gf ) = 



ab 

where we have also introduced their diagrammatic 
representation. 

In the following, for simplicity we will work in 
the symmetric phase where the macroscopic field 
expectation value vanishes, i.e. (po(k) — 0. As a 
consequence, also all three-vertices vanish for the 
considered theory with interaction (l45t and using 
our above example we get the comparably compact 
expression: 

^Gf^l 



c(x)6<f> d (y) 



_ /^R r-</><P<P</> /^R 
~ ^k.ae 1 k,efcd U k.fb 



0,0 = 



, :p r## /-R 
+ lf k ,ag^ k ,g hcd ^ k ,hb • 



We will use a diagrammatic representations, 
where the retarded and advanced cutoff functions 
R^' are denoted by the insertion of a cross, i.e. 

b 



K,ab 



a 



R 



k,ab 



a - — x - - b 



and the proper vertices, indicated by the full dot, 
are given by 



k.abcd 



X,- 



k.abcd 



v 
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and equivalently for the remaining four-vertices. 

For the exact flow equation for the effective av- 
erage action the one-loop form is written diagram- 
matically as 



' ;:(*>• O 



Similarly, using a compact notation for the deriva- 
tives the flow equation for the two-point function 
takes the form 

a b a b 



0-0 



a b a b 

and for the four-point function we have 

p(4) 

k,abcd 



abed abed 


+ A-A + a4\ 

abed abed 


+ At- A + fa^k 

abed abed 

+ h~f\ + hrK 


abed abed 

+ /TTv + fork 


abed abed 
+ P(a, b, c, d) | 

H + + 

abed abed 


abed abed 
abed abed 



Here, P(a, b, c, d) denotes all possible permutations 
of the indices on legs of the respective diagrams. 
The diagrammatic representation of the flow equa- 
tions for n-point functions follows the standard 
construction rules: draw all combinations of prop- 
agators Gf' A and Fj, and vertices Tr' and attach 
the appropriate symmetry factors, taking into ac- 
count that certain diagrams may either vanish or 
are identical. We emphasize, however, again that 
in the diagrammatic representation the propagators 



are evaluated at an extremum of the effective aver- 
age action. 

We may simplify the construction rules of the 
flow equations even further by introducing the 
derivative operator d k , 



d k T k [4>,4>] =Tr 



Rf—z+Rf — r 

k 6Rf k 6Rf 



(65) 



With this we may reduce expressions like 

(rf + Rff 1 Rf (rf + Rff 1 = -~d k (i* + Rf)' 1 , 

to a simple form. For instance, the flow equation 
for the two-point function can be written as 



K2) 

k,ab 



;-iQ . Q • Q 



a b 



Equivalently, for the four-point function H we 
have 



K4) 
k,abcd 



-*m y:y. 



+ 

d b 



a 

+ 

d b 



b' *S V 

P{a, b, c, d) 



fl V~V c fl v~V £ 



4M 





abed abed abed 

5. Solving truncated flow equations 

By calculating the functional derivatives of the 
effective average action, we extracted above the 
flow equations for «-point functions. We will use 
these flow equations in the following to study the 
behavior in the vicinity of possible nonthermal 
fixed points. Fixed points correspond to trans- 
lationally invariant scaling solutions for n-point 
functions both in space and time. Therefore, we 
are interested in the limit to —> -oo for which the 
dependence on the details about the initial condi- 
tions encoded in g(to) are lost. As discussed in 
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section [3] conserved quantities will play an im- 
portant role and we will not impose a fluctuation- 
dissipation relation to be able to describe nonther- 
mal fixed points. 

5.1. Stationarity condition 

The presence of nonthermal scaling solutions 
was discussed in section [3] based on a stationar- 
ity condition for the nonequilibrium time evolution 
equations. An equivalent (scale-dependent) sta- 
tionarity condition can be obtained from the func- 
tional renormalization group, where it appears as a 
nontrivial identity relating the various second func- 
tional derivatives of Fa[0, </>]. For spacetime trans- 
lation invariant systems, it is very convenient to 
consider the correlation functions in Fourier space. 
We start by writing down the following identity in 
momentum space, 

= -iGl ac (p)lf cd G£ de (p) 

x {( G *)>H G *)>)} " (66) 

In order to interpret further the different combina- 
tions of terms appearing in this equation, we write 

(2) 

the two-point functions I"V in the form 



k I 



o -(<?) 

iD _1 \ / 
iD 1 



^1 



where iD 1 is the free inverse propagator. The sec- 
ond line defines the self-energies 2^ ' A ' ; 
turn are given by 



which in 



-'k.db 



k,ab 



^k.ab ~ ^k,ab ^k.ab ~ 1 tab 1 k,ab 



employed in section [3] it essentialljQ states that 
'gain terms' equal 'loss terms' for which station- 
arity is achieved I22I1 . This aspect will be analyzed 
in detail in section [6] 

Of course, the condition d67l i is trivially fulfilled 
if the fluctuation-dissipation relation holds, i.e. in 
thermal equilibrium where 



Ft"\p)=-i[n BE (p°) + ^yr\p), 

^ m (p) = -i{n B E(p°)+%f ci> (p), 



for the propagators and self-energies I122II . These 
relations will not be assumed in the following. 

Using the representation of the self-energies in 
terms of two-point functions we may immediately 
write down the flow equations for the statistical 
component Sf , which is given by 

= \~ dk { Q + Q + Q 

a b a b a b 

reading in momentum space 



The spectral self-energy satisfies the flow equation 



k,ab k,ab k,ab 



= -hdi 



{ Q ■ C Q 

a b a b a b 

-Q-G-Q 

a b a b a b 



where we have used that Rf = Rf for the consid- 
ered class of cutoff functions ( f59b - With these iden- 
tifications, the identity d66l l can be written in terms 
of the self-energies as 

Kab^PkMiP) ~ Fk.ab(p) Z P kM (p) = . (67) 

This equation is well-known in nonequilibrium 
physics. In the language of Boltzmann dynamics 



4 Using the Wigner coordinates employed in section fj] the 
nonequilibrium time evolution of the statistical function is given 



by [2 

2p"<V F ab (X, p) = i [l. F ac p cb - F ac lT b ) (X, p). 

The condition (67) is related to this by taking the trace at a sta- 
tionary point, where the correlation functions become indepen- 
dent of X. 
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5.2. Resummed 1 jN -expansion 

In general, flow equations for n-point func- 
tions are expressed in terms of (n + 2)-point func- 
tions, and to find approximate solutions one has 
to truncate the infinite hierarchy of coupled flow 
equations. Here, we use an ansatz for the four- 
point functions, which is equivalent to a resummed 
large-Af expansion of the 2PI effective action to 
next-to-leading order (NLO) lIHoll . 

Using 0{N) symmetry, we can always write the 
two-point functions as, for instance, 

rf>,>0 = rf(*,)0^. 

The four-point functions can be decomposed into 
the contributions from the different channels, e.g. 
forrf # , we havelH 



to*.***) 

= if ' *\x,z)5 (d+l \x - y)6< d+l Xz - w)6 ab 6 cd 
+ lf^(z,y)6 (d+l \z - x)6 (d+l \y - w)6 ac S bd 
+ I?*(y, z)6 (d+l) (y - z)6 {d+l \x - w)6 bc 5 ad 



In the resummed large-A^ expansion we write these 
individual contributions diagrammatically as 



n pm {x,y)= x) 



where the full blobs denote the resummed bubble 
chain with full propagators G R,A , F, F inserted on 
the internal lines connected by the bare vertices 
given by 



A 

'3N 



A 
12N 



For instance, 



>•< = x - m >Ck 
= x - ^ >*OK 

Each vertex contributes a factor 1/Af and each 
closed loop gives a factor of N = 5 a b5b a - Thus, 



all shown diagrams contribute at the same order, 
and there are no other diagrams that contribute at 
this order in the expansion. Notably, all n-point 
functions rl"' with n > 4 are of higher order in 
the large-A^ expansion and the infinite hierarchy of 
flow equations is closed at the level of the four- 
point functions 11171 Il8ll . Diagrammatically this 



can be expressed in terms of the flow equation 



>< 



iN d k 



which has only the four- vertex Yv' appearing on 
the r.h.s. We also note that the scale derivative 
represents a total derivative in this approximation 
Ill7l,ll8l 16111 . which can directly be understood from 
our derivation of the exact flow equation with the 
help of the 2PI effective action in section [4] 

Similarly, one obtains for the other four-vertices 



iN 



+ ( - iN)2 ;OC< 



- iN 



>(X 



and 



t > ( ^{ X)-.~X ■ >(X)< 
+ XXX + XXX 
XXX XX.X 
4"{ MX + XX 

+ xx 



The loop contributions appearing in the chain of 
bubbles diagrams are given by the expressions 



F k (p-q)G,' (q) 



(68) 
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and 



(p) = ^ J [F k (p - q)F k {q) - \pk{p ~ q)Pk(q)] ■ 

(69) 

In terms of these one-loop expressions the four- 
point functions read 



and the spectral component 



y m ( = rf*{i-nf + (nf) 2 + .. 

_J_ d££ n A 
3N 3N k ' 



A 

'3N 



-nf 



3N 



K) 2 



+ . . 



A 
V2N 



+ UN k 



where we have defined the momentum-dependent 
effective coupling 

,1 



'WO) = 



(70) 



[i + n>)][i + n>)] ' 

to express the set of resummed diagrams. In or- 
der to observe the equivalence with the previous 
expressions above, we note that, e.g., for the set 
of diagrams contributing to the vertices and 
F^^, we have 
1 1 



+ . . 



i+n»i + n» 

= l - nf(p) - Ufa) + (nf(p)f + (n£( P )f 

where mixed products of the type Tlf(p)Ylf(p) van- 
ish in the expansion, since there are no vertices 
r** that would allow for that particular combina- 
tion of retarded and advanced one-loop diagrams 

n R,A 

Plugging the vertices into the flow equation for 
the self-energies Jf k ' F and integrating the total scale 
derivative gives the final result for the statistical 
component 



- i n k(P ~ i)Pk(q) 



Writing these expression explicitly in momentum 
space, we have 

K (p) = ~ f o Aeff - k(p ~ q) K (p - 4) F *te) 

(71a) 

z p k (p) = -jjj J A e «,k(p - q) { n * G> - q)Pk(q) 

+ Yf k (p-q)F k (p)} , (71b) 

which has the structure of a two-loop self-energy, 
however, with a momentum-dependent coupling 

^eff.A- 

To summarize, we have found that in the large- 
N expansion to NLO the infinite hierarchy of flow 
equations is closed on the level of four-point dia- 
grams. These can be solved directly by integrating 
the total scale derivative d k . Thereby we obtain the 
full expressions for the self-energies and vertices 
to this order. In the following section we will in- 
vestigate the scaling behavior of the self-energies 
that enter the stationarity condition (l67t . That way 
we can extract the scaling exponent k characteriz- 
ing different types of fixed points. 

6. Strong vs. weak wave turbulence 

6.1. Nonperturbative stationary transport 



From the integrated self-energies (171 at and 
( 171bl ) we can directly classify the scaling solutions 
in the limit k —* 0, where the regulator is sent to 
zero. In order to better compare the nonperturba- 
tive aspects of these results to the perturbative dis- 
cussion of section[3] we write without loss of gen- 
erality for k = 0: 



F{p) = -iL{ P )+ l -y 



(72) 



and the presentation follows to a large extent 
Ref. ||44|]. Equivalently to ( fT8b . the function n{p) 
depends on the four-momentum p = (p°,p), in 
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contrast to the case of thermal equilibrium (see sec- 
tionQ}. It satisfies the symmetry property 

n(-p) = -(n(p) + l) , 

which follows from F(—p) = F(p), and p(—p) - 
-Pip). 

We then write the stationarity condition (|67| | in 
terms of n(p) and the spectral function p(p) using 
the identity 

If{p)F{p) - •L F {p)p{p) 

= j(^(p)-^(p))(F(p) + ip(p)) 

- i{lF{p) + ^E p (p))[F(p) - ip(p)) , (73) 
where by (l72l i we have 

F(P) + ^P(P) = ~in(p)p(p) , 
F{p)- l -p{p) = -i(n(p)+\)p{p) . 



From the self-energies (17 lab and (I71bt we con- 
struct the linear combinations 

Yf(p) + l -If(p) 

I AeffiP — k{) 

18Njl,2 

x (F(p -h- k 2 ) + l -p(p -ki- fc 2 )) 

X ( Fl 4 Pl )(^ 2 4 P2 )' 
that enter the stationarity equation (l73~t . In terms 
of n{p) this reads 

x (n(p -ki- k 2 ) + I) p(p -k { - k 2 ) 
x (n, + l)pi (n 2 + l)p 2 , 
S F (p) + l -V{p) = -» J A^(p - fc) 

x n(p - ki - k 2 )p{p - ki - k 2 ) 
x «ipi n 2 p 2 . 

Finally, putting everything together we obtain the 
following form for the stationarity condition: 

-i(EPF-?Fp){p) 

= -4X 2 ,3 (2 ^ +1)( ^- &1 - fc2 - &2) 
X ,leff(P ~ ("1 + 1)(« 2 + 1)("3 + 1)",, 

-nxnjn-iiyip + l)|p I p 2 p 3 Pp . 

(74) 



We may bring this expression into a form which 
can be directly compared to kinetic or Boltzmann 
descriptions by mapping onto positive frequencies 
p°. In particular, we consider the 'collision inte- 
gral' 

^(S'T-E^)(p°,p) S C NL >), 

which will allow us to relate this discussion to 
the presentation in section [3] Here, the upper in- 
dex NLO indicates that this is accurate to next-to- 
leading order in the large- N expansion, including 
processes to all orders in the coupling constant, in 
contrast to the perturbative discussion in section [3] 
After performing the positive frequency integral, 
we get 

C NL0 (p) = 

JdF 2 ^ 2 + 1) («i + 1) n 2 n 3 

- n p n\ (« 2 + 1)(«3 + 1) j 
+ J" difl 3 {(n p + 1) («, + 1) (« 2 + 1) «s 

- n p nin 2 (n 3 + 1) J 
+ JdT® 3 [{n p + l)n 1«2«3 

-« ; ,(«1 + l)(n 2 + D(«3 + 1)} 
+ JdT »4 \{n p + 1) (m + 1) (« 2 + 1) («3 + 1) 



-n p ni«2«3 



(75) 



with the 2 3 = 8 contributions from the different 
orthants in frequency space and 



= 71^ f (2n) d+1 * d+l) (P + ki-k 2 -k 3 ) 

LOl\ J 1,2,3 

x [Aesip + h) + A eB (p - k 2 ) + Aesip - k 3 )] 

X PlP2P3Pp , 



(a) 
l<->3 



= -i- f (2n) d+1 ^ +l \p + h+k 2 -k 3 ) 

J 1,2,3 

x [A eB (p + h) + A eS (p + k 2 ) + Aesip - k 3 )] 

x plp 2 p3Pp , 
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.(b) 
1h3 



-i- f (2n) d+i 6^ l \p - kl -k 2 -k 3 ) 

ioW J 1,2,3 

X A eff (p - ki) p 1 p 2 p3Pp , 



= 77^7 f (> W +1 ^ +1) (/> + *i + * 3 ) 

J 1,2,3 

X ^eff(P + h)p\P2PiPp ■ 



Here, the (>) sign on the integrals indicates, that 
the integrals run over positive frequencies, i.e. 



r(>) r a 
J 1,2,3 Jo 



2 * 11 2 « 



]/ 


r-r d% 




Ui.3 ( 2 -) d J 



From the Boltzmann description of the weakly 
interacting theory we clearly recognize from the 
first contribution above 2 <-» 2 scattering processes, 
however, with an effective coupling A^(p). This 
momentum-dependent coupling is a consequence 
of the summation of an infinite number of pro- 
cesses, which will be crucial in order to be able to 
discuss the nonperturbative regime of strong tur- 
bulence at low momenta as is explained in the 
following. All other processes are 'off-shell' and 
turn out not to play an important role in this con- 
text HE! H. 



6.2. Scaling solutions 

We want to study the nonequilibrium steady 
states and their scaling properties. Along the lines 
of our discussion in section[3] we can extract these 
properties from a scaling analysis. For a discus- 
sion in terms of a direct determination of the prin- 
cipal zeros of the collision integral (175) . we refer 

toRefs. HE HI! 

From the quantities entering the collision inte- 
gral d75l l it is clear that we need the scaling proper- 
ties of n(p), the effective coupling A e s(p), and the 
measure. For that purpose, we consider first the 
scaling behavior of n(p). Extending the discussion 
of section [3] we also take into account a possible 
dynamic critical exponent z different from one and 



5 Ref. fT3l tacitly assumes the absence of particle number 
changing processes for the derivation of the part icle cascade. 
For a proper discussion of this aspect see Ref. Iia. 



a nontrivial anomalous dimension rj following sec- 
tion [2] With the scaling ansatz given by ( 111a) and 
( II lbb for the statistical and spectral function, re- 
spectively, we have 

n(p°,p) = s K+,1 n(s z p , sp) , 

again assuming that n(p) » 1 /2. 

For the scaling analysis of the momentum- 
dependent effective coupling A s g(p), we use 

(n A ( P )f = n A (- P ) = n R ( P ) , 

to write d70t as 

,1 



U + n R (p)p 



(76) 



To extract its scaling, we need the scaling behavior 
of the one-loop diagram 

n R (p)= j f F(p - q}G K (q) , 

which is obtained from the statistical propaga- 
tor F and the representation of retarded propaga- 
tor in terms of the spectral function G R (x,y) = 
p(x,y)0(x° -y°), i.e. 

G\p ,p) = s 2 -»GXrp°,sp). 

Thus, we have 



n R (p) = | jF(p - q)G R (q) 



f s 2+ «F(s z (p -q°),s(p-q)) 

Jq 



x s^G^sY, sq) , 

for the one-loop diagram n R . By an appropriate 
rescaling of the measure of this one-loop diagram, 
i.e., taking q° — > s~ z q° and q — > s l q, it can be 
brought to the form 

n R ( P ) = s A n R (s z p°,sp), (77) 

with the scaling exponent 

A=4 + K-ij-z-d. 

For positive A > one observes that Tl R (p) » 1 
for sufficiently low momenta, such that fluctuations 
become important in the infrared. Inserting this re- 
sult into the expression d7o*l i for the effective cou- 
pling, we finally get 

A > : A eS (p°,p) = s- 2A A eff (s : p°, sp) . 
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In contrast, for the case A < the scaling of the 
effective coupling is trivial. This case will be rele- 
vant at sufficiently high momenta, where Il R (p) is 
small. Then, from (|76*| | it follows that the effective 
coupling is essentially equivalent to the perturba- 
tive coupling, 

A < : A eS (p) - /I • 

It remains to determine the scaling behavior of 
the measure J dT. We consider, for instance, the 
2 <-> 2 processes: 



dr 



2<->2 



_j_ rfe! n ^ 

18N Jo ( In . = \ 1 3 2n 



(27T) d+1 6« +l \p + k 1 -k2-k 3 ) 



x [A e g(p + + A eS (p - k 2 ) + Aeftip - k 3 )] 

*PlP2P3Pp ■ 

From the scaling analysis, we obtain 
J dr2n2(s z p°, sp, s z k®, ski, s : k®, sko, s : k®, ski) 



where the momentum dependencies in the measure 
are indicated explicitly. With a positive exponent 
for the one-loop diagram (l77t . A > 0, the scaling 



of the measure can be written as 



-2k-2t]-z 



In 



contrast, for A < 0, where the coupling is given 
by /l e ff — A, one finds the scaling ~ s - i :- 2it+& - 4 i ^ 
The same scaling properties are also obtained for 
the remaining measures for the 'off-shell' 1 <-> 3 
and <-> 4 processes, which may be neglected for 
the following discussion. 

Putting everything together, in the infrared 
where strong fluctuations dominate the dynamics, 
we obtain 



A>0: 



C NLO (p) 



K+rj-z fNLQ 



C"^{sp) . 



In contrast, in the case where the scaling of the cou- 
pling is trivial, i.e. /t e ff(p) — A, relevant at high mo- 
menta we have the scaling behavior 

A < : C NL0 (p) = s 3 "-i- 3: - 2d+ *C NLO (sp) . 

First, we consider the implications of the scal- 
ing behavior in the high momentum region where 



A < such that A e fi(p) A. In this case we re- 
cover the perturbative discussion of section [3] and 
the momentum integral 



d\p\\p\ d - l C Nh0 (p) 



describes the particle flux through a sphere of ra- 
dius k. Performing the scaling transformation, we 
obtain 



A<0: k = d + z + 



for the scaling exponent k. Setting the anomalous 
dimension 77 = and the dynamical critical expo- 
nent z = 1 one may immediately verify that this 
reproduces the same weak wave turbulence expo- 
nent k — I for the particle cascade in d — 3 spatial 
dimensions as in section[3] 

In contrast, in the infrared scaling region with 
positive values of A, we obtain the exponent 

A>0: K-d + z-r], 

for the particle flux. This scaling behavior in the 
nonlinear, low-momentum regime is also known 
as strong turbulence, which is displayed schemati- 
cally in Fig. Q] in section [3] Of course, we can also 
extract the scaling behavior associated with the en- 
ergy cascade. The scaling solutions for the rela- 
tivistic scalar A^-component theory, where z = 1 
and 77 = are summarized in the table. The pre- 
dicted values for the exponents have been tested 
using classical-statistical simulations on the lattice 
in various dimensions J5I S ft Similar studies 
have also been performed for nonrelativistic scalar 
or in the context of non- 



theories (l2j \l% |45[ [62 



abelian gauge theories 



7. Summary and outlook 

In these lectures we have discussed different as- 
pects of the nonequilibrium functional renormal- 
ization group. As an example, we considered the 
physics of wave turbulence in relativistic scalar Af- 
component field theory. We have seen that the stan- 
dard treatment based on Boltzmann equations can 
only describe the perturbative regime of the associ- 
ated stationary transport of conserved charges. In 
contrast, the functional renormalization group can 
efficiently describe the perturbative regime as well 
as the nonperturbative physics at low momenta. 
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Particle cascade 


Energy cascade 




Strong turbulence 


k = d + 1 


K = d + 2 




Weak turbulence 


k — d - 4 


k — d — | 


quartic interaction 




k — d - 2 


K = d-\ 


cubic interaction 



We have approximately solved the renormaliza- 
tion group equations using a (2PI) resummed large- 
N expansion to NLO. This approximation is par- 
ticularly suitable, since it provides an accurate de- 
scription of the physics at not too small N, while it 
is still analytically treatable. In particular, the scale 
derivatives are total derivatives in this approxima- 
tion, which can be trivially integrated. This al- 
lowed us to efficiently discuss the close relations 
and differences of the renormalization group ap- 
proach with kinetic theory and 2PI effective action 
techniques in nonequilibrium physics. 

The renormalization group provides powerful 
nonperturbative approximation schemes beyond a 
resummed large-Af expansion, such as derivative 
expansions or expansions in powers of fields with 
momentum-dependent vertices which we have not 
considered in these lectures. While this has been 
explored to a large extent mainly in Euclidean 



spacetime for physics in thermal equilibrium J32 
4211 . much less is known in the context of far 



from-equilibrium problems in quantum field the- 
ory. Here the functional renormalization group 
may give important insights also for the nonequi- 
librium dynamics of, in particular, non-abelian 
gauge theories, where suitable large- N techniques 
are difficult to implement. 
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